GENERALIZING DODGSON'S METHOD: A "DOUBLE-CROSSING" APPROACH TO 

COMPUTING DETERMINANTS 



Abstract. Dodgson's method of computing determinants was recently revisited in a paper that appeared in 
the College Math Journal. The method is attractive, but fails if an interior entry of an intermediate matrix has 
the value zero. This paper reviews the structure of Dodgson's method and introduces a generalization, called a 
"double-crossing" method, that provides a workaround to the failure for many interesting cases. 



1. Introduction 

Algebra students learn this simple pattern for the determinant of a 2 x 2 matrix: 



X 



d 



ad — be. 



A similar pattern exists for 3 x 3 matrices, but to compute the determinant of larger matrices, the student 
must learn something a little more complicated. Most students learn to compute determinants by expansion 
of minors, first developed by Laplace. Some students also learn to compute determinants by triangularizing 
the matrix. Both methods are effective, and triangularization is efficient, but 

• hand computations frequently lead to many mistakes; 

• expansion of minors is tedious; and 

• triangularization can turn a matrix of integers into a matrix of fractions. 

In 1866, the Rev. Charles Lutwidge Dodgson 1 , developed a conceptually simple method to compute deter- 
minants [2]. Dodgson's method iterates the familiar 2x2 formula; although it uses division, matrices with 
integer entries do not turn into matrices with rational entries. 

Example 1. Given the matrix 

10 1 
1 3 1 
1 1 

we set A^ — A. To compute A^ 2 \ compute the determinants of the 2x2 contiguous submatrices of A^: 



A (2) 



/ 


1 




1 


\ 




1 3 




3 1 






1 3 




3 1 




V 


1 




1 1 


/ 



To compute A^\ repeat the computation for A^ 2 \ but divide by the center entry of A^: 



V 



3 -3 
1 2 

3 



(3). 



Now go back and compute det A using your favorite method (pattern, expansion of cofactors, triangular- 
ization, etc.). You will find that det A = 3. 0> 

Dodgson's method is quick and conceptually simple. In general, we can describe the method in the follow- 
ing way: 

• Let A(") be the n x n matrix given. 



Dodgson is better known to children of all ages as Lewis Carroll, author of "Alice in Wonderland" and "labberwocky" 

1 
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• For each k — n — 1, n — 2, . . . , 1: 

- Let B^ be the k x k matrix of determinants of contiguous 2x2 submatrices of A^ k+1 \ 

- If£: = n-l,letA« =B(*). 

- If A: < n — 2, let be the k x k matrix whose (z',/)-th element is the (i,j)-th element of B^ 
divided by the (z',/)-th element of the interior of A^ k+2 \ (The interior of an n x n matrix M is 
the (n — 1) x (n — 1) submatrix whose (z',/)-th element is the (z + l,j + l)-th element of M). 

• The singleton element of is the determinant of A. 

Dodgon's method is an example of a condensation method to compute determinants; each iterate A^> is 
a condensation of the previous iterate A^ k+1 h Other condensation methods appear in [1]. If Dodgson's 
method terminates successfully, it computes the determinant of an n x n matrix using 



4 (n-l) 2 + (n-2) 2 + 



+ v 



i-n-n 



multiplications, subtractions, and divisions, or O (n 3 ) operations in Q. This is not bad, especially consid- 
ering that all the divisions are exact, so no fractions are introduced; Bareiss' algorithm, a better-known 
fraction-free method of computing determinants, also performs O (n 3 ) operations in Q. [3, 4] 

There's the rub, though: division presents Dodgson's method with a huge drawback. (It also presents 
an obstacle with the usage of Bareiss' algorithm.) What's so bad about division? 

Example 2. Swap the first two rows of the matrix of Example 1 to obtain 



M 



We know from the properties of determinants that det M = — det A — 
What happens when we use Dodgson's method? First we compute 




-3 3 
1 -1 



To compute we divide the determinant of M^ 2 ' by the interior element of . But the interior element 
of M^ 3 ) z's zero! <> 

In general, Dodgson's method fails to compute the determinant of a matrix A whenever a zero appears 
in the interior of A W for any k > 3. This can happen even if no zeroes appear in the interior of A. 

A workaround discussed in [2] swaps rows of the original matrix in such a way that zeroes are moved out 
of the interior. For example, if you swap the top two rows of M in Example 2, you return to A of example 1, 
for which Dodgson's method worked fine. However, there are two drawbacks to this workaround. First, 
swapping rows may well introduce other zeroes into the matrix, and it isn't easy to predict this from the 
outset. Second, swapping rows simply won't work for some matrices. 

Example 3. No combination of row or column swaps will allow Dodgson's method to compute the deter- 
minant of 

/ 1 3 \ 
0-101 
1 12 
\ 2 1 j 

because there will always be a zero in the interior of N. C> 

Will a different workaround of Dodgson's method work for N? Yes! We describe such a method in 
Section 2. The reader will see quickly why we call it a "double-crossing" method. 



N 



2. A "DOUBLE-CROSSING" METHOD 



Before describing the method, we illustrate it using the matrix from Example 2. 
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Example 4. Recall 



(2.1) 




and M (2) = 



-3 3 
1 -1 



The zero in the interior of causes Dodgson's method to fail when computing My> . 



(3) 

Above that zero is a non-zero element, M\ 2 = 3. We will divide by this element instead, but this requires 

us to re-compute in a slightly different manner. Cross out the first row and second column of M^ 3 ' 
(the ones containing 3). We are left with the 2x2 complementary matrix 



M* 



Recall that we compute the minor of an element of a matrix by (again) crossing out the row and column 
containing that element, then taking the determinant of the remaining submatrix. Consider the matrix M' 
of minor of elements of that correspond to the elements of M*; that is, 





/ 


3 1 




1 3 


\ 






1 1 




1 




M' = 
















3 1 




1 3 






V 


1 




1 


/ 


M( 2 ) = 


M' = 


[2 P 
13-3, 





Put 

(2.2) 

conclude by dividing detM' 2 ^ by the non-zero element of M^ 3 ' that we identified earlier: 

M( 2 ) 

3 



-6-3 
3 



(-3). 



As noted in example 2, det M = —3. 

Notice, by the way, that after a row and column swap all but one of the values in M^ 2 ' in line (2.2) are 
the values of in line (2.1). We will say more about this later. 0> 

The generalization we have presented of Dodgson's method preserves its "spirit", inasmuch as we com- 
puted all determinants by condensing 2x2 contiguous submatrices. The example shows why we call it the 
"double-crossing" method: 

• we found a non-zero element adjacent to the zero element; 

• we crossed out its row and column, obtaining the complementary matrix M*; 

• for each element in M* , we compute its minor by 

- (again) crossing out the row and column of its location in and 

- using the determinant of the remaining matrix to compute the element of . 

When a zero appears in an intermediate matrix AW (where n > k > 3), we cross out rows and columns 
that make up a submatrix of A^ n \ Let's look at a matrix where a zero does not appear in the interior of A, 
but does appear in the interior of an intermediate matrix. 



Example 5. Let 



A = 



( 1 1 1 \ 
5 3 1 
13 2 11 
1110 
V 2 2 1 j 
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Using Dodgson's method, we compute A^ 5 \ and then 

A( 3 ) = 



-4 
-2 
4 



-2 

-4 



We encounter a zero in the interior! Above it is a non-zero element, —2. Notice that it lies in the top row 
and central column of 

To compute A^ 2 \ do the following: 

• cross out the top three rows and central three columns of A^ (the original matrix); 

• this gives us a 2 x 2 complementary matrix M*; 

• for each element in M* , 

- cross out the row and column of its location in A^\ 

- compute the resulting minor, and 

- put it into the corresponding location in A^ . 
Following these instructions, we have 



and 



A (2) 



/ 





1 





1 




1 





1 





\ 




5 


3 


1 










5 


3 


1 






3 


2 


1 


1 




1 


3 


2 


1 









2 





1 




2 





2 












1 





1 




1 





1 









5 


3 


1 










5 


3 


1 






3 


2 


1 


1 




1 


3 


2 


1 




V 


1 


1 


1 










1 


1 


1 


/ 



We compute the determinants for A^ using Dodgson's method — and again, the bottom two values turn 
out to have values that you would obtain in the two top rows by ordinary condensation of A^\ while the 
others are different: 

^4-4 

We now compute A^' by dividing the determinant of A^ by the non-zero entry of A^ that we identified 
above: -2. We obtain 



= (4), 



and in fact det A — 4. O 



So far we've been using a special case of the double-crossing method. Theorem 6 describes this special 
case precisely. 

Theorem 6 (Double-crossing method, special case). Let Abe ann x n matrix. Suppose that we try to evaluate 
det A using Dodgson's method, but we encounter a zero in the interior of A^ k \ say in row i and column j of A^ k \ If 
the element a. in row i — 1 and column j of A^ is non-zero, then we compute A^ k ~^and A( fc ~ 2 ) as usual, with the 
following exception for the element in row i — 1 and column j — 1 of A^ k ~ 2 >: 

• let £ = n — k; 

• identify the (£ + 3) x (£ + 3) submatrix A whose upper left corner is the element in row i — 1 and column 
j - 1 ofA^; 

• identify the 2 x 2 complementary matrix M* by crossing out the (1 + 1) x (l + 1) submatrix of A whose 
upper left corner is the element in row 1 and column 2 of A; 

• compute the matrix M' of determinants of minors of elements ofM* in A; 

• compute the element in row i — 1 and column j — 1 of A^ k ~ 2 > by dividing the determinant of M 1 by a. 
We can use Dodgson's method to compute the intermediate determinants ofM'. O 
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A proof of correctness appears in Section 4, after we explain why Dodgson's original method works 
correctly. 

One can generalize Theorem 6 so that the non-zero element appears immediately above, below, left, right, 
or catty-corner to the zero; that is, the non-zero element is adjacent to the zero: see Theorem 13 on page 13. 
If the zero appears in a 3 x 3 block of zeroes, then the double-crossing method will not repair Dodgson's 
method, although additional strategies may be possible. 

Before proceeding to the next section, we encourage the reader to go back and examine how Theorem 6 
describes what we did in the examples of this section. We conclude by applying the double-crossing method 
to compute the determinant N from Example 3. 

Example 7. Recall from Example 3 



N = 



( 1 



1 

V o 





-1 
1 

2 



\ 
1 



1 J 



Put N( 4 ) = N: we have one zero element in the interior, at = (2,3). As described in Theorem 6, this 
corresponds to element (1,2) of N^ 2 '. The other three elements of can be computed as usual: 

/ -1 3 3 \ 
= 1-2-2 — > 
\ 2 -4 2 ) 

We use the double-crossing method to compute the other element of N^ 2 h The interior zero appears in the 
original matrix, so t — 0. We choose the 3x3 submatrix whose upper left corner is the element in row 
i - 1 = 1, column; - 1 = 2 of 

/ 3 
A=\ -1 1 
\ 1 2 

Cross out the lxl submatrix in row 1, column 2 of A and identify the complementary matrix 

-1 1 
1 



M* = 



Compute the corresponding matrix of determinants of minors in A 



M' 



( 



V 



/ 



The determinant of this matrix is 9; after dividing by the nonzero N^j = 3 we have 



NO 



We can now conclude by computing 



(3). 



In fact, the determinant of N is 3. Q> 

The most computationally intensive part of the double-crossing method is that of computing the ma- 
trix of determinants of minors, but many of these have already been computed. It likely is not clear to the 
reader at present how to identify these, but as we explore the mechanics of Dodgson's method and the 
double-crossing method, we will see that we can predict exactly which parts of the matrix of cofactors need 
recomputing, and which can be copied from previous work. 
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3. Why does Dodgson's method work? 

Just as Bareiss' algorithm relies on a well-known theorem of Sylvester, Dodgson's method relies on the 
following theorem of Jacobi. [2] 

Theorem (Jacobi's Theorem). Let 

• Abe an n x n matrix; 

• Manmxm minor of A, where m < n, chosen from rows i\, ii, ■ ■ ■ ,i m and columns ■ ■ ■ ,j m , 

• M' the corresponding m x m minor of A', the matrix of cof actors of A, and 

• M* the complementary (n — m) x (n — m) minor of A. 

Then 

detM' = (detA)™' 1 ■ detM* ■ (-l)^ if+it .O 
Section 4 gives a proof of the new method. 

We adopt the following notation. If A € K mx " denotes a matrix, then A' denotes its matrix of cof actors. 
Both Dodgson's method and the double-crossing method make use of submatrices of a matrix, which for a 
given matrix A we denote in the following fashion: 

Ai.,.j / k...£ rows i, i + l,...,j and columns k, k + 1, . . . ,1 of A; 
A',)c...e row i an d columns k, k + 1, . . . , I of A; 
A'.../,jt rows i t i + 1, . . . and column k of A. 

Example 8. Recall from Example 2 

/ 1 3 1 \ 
A = 1 1 . 
\ 1 1 / 

Its matrix of cofactors is 



A' = 



f 





1 




1 


1 




1 





\ 




1 


1 







1 







1 






3 


1 




1 


1 




1 


3 






1 


1 







1 







1 






3 


1 




1 


1 




1 


3 




\ 





1 







1 




1 





/ 



while 

A!... 2/ i...2 = ( 1 o ) and A3 ' 2 - 3 = ( 1 1 ) • ° 



To illustrate the relationship between Dodgson's method and Jacobi's Theorem, consider a generic 4x4 
matrix. 
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FIGURE 3.1. Diagram illustrating how Jacobi's Theorem is applied to Dodgson's method, 
given the 4x4 matrix of Example 9. 











^2..3,2..3 











Consider the generic 4x4 matrix from Example 9. The final condensation from Dodgson's method 
gave us 

^(1) _ / 1^1. .. 3,1.-3 1 1^2... 4,2.. .4 1 ~ |-^2...4 / 1...3l|-^1...3 ) 2...4| 
V \ A 2... 3,2.. .3 1 

For Jacobi's Theorem, select the 2x2 minor of the corners of A 



M 



M,i Ma 

A-4,1 ^4,4 



Its complementary minor is M* = (A2...3,2...3). The matrix of cofactors is 



M' 



A l,l ~ A 'l,4 



-A' A' 
/i 4,l ^4,4 



where A'^ = \A 2 ..a,2..a\, A' u = \A X „ 3A „ 3 \, A' u = \A 2 .,a,i...3\, and A^ = |A L .. 3 , 2 ...4|- Byjacobi'. 
Theorem, 

detM' = (detA) 2 " 1 - detM* ■ (-l) 1+1+4+4 
detM' 



detA 
detA. 



detM* 

|^1...3,1...3| |^2...4,2...4| — |^2...4,1...3l |^1...3,2...4| 
| A2.. .3,2.. .3 1 

This is precisely the final step in Dodgson's method as applied to A. Notice that the negatives in M' 
cancel when computing the determinant. 



Example 9. Let A be a generic 4x4 matrix; we show how Dodgson's method applies Jacobi's Theorem. 
The first two condensations by Dodgson's method produce 



(3.1) 



A (3) 



I l^l...2,1...2| |^1...2,2...3| l^l...2,3...4| ' 

IA2...3/I...2I |A2...3,2...3l |A2...3,3...4| 
\ |^3...4,1...2| IA3... 4,2.. .3 1 IA3...43...4I ) 



and 



1...24...2 



|Ai...2,2... 3 | 



|^2...3,1...2| IA2... 3,2.. .3 1 



I Ai.. .2,2.. .3 1 |Ai...2,3...4| 
1^2.. .3,2.. .3 1 |^2...3,3...4| 



"2,2 



3,1...2| 



2...3,2...3| 



|^3...4,1...2| IA3.. .4,2.. .3 



|^2...3,2...3| 
l^3...4,2...3l 



«2,3 



l^2...3,3...4| 
l^3...4,3...4| 
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To see how corresponds to Jacobi's Theorem, consider the upper left 3x3 submatrix of 

(01,1 «1,2 fl l,3 

«2,1 fl 2,2 «2,3 

fl 3,l fl 3,2 fl 3,3 

Cross out row 2 and column 2 of A, obtaining 

M—( Su fll ' 3 
V fl 3,l fl 3,3 

its complement in A is M* — (#2,2)- The 2x2 submatrix of A' corresponding to the cof actors of M in A is 



M' 



/ 


«2,2 #2,3 




«2,1 «2,2 


\ 




fl 3,2 fl 3,3 




fl 3,l fl 3,2 






«1,2 «1,3 




«1,1 «1,2 




V " 


«2,2 «2,3 




fl 2,l fl 2,2 


/ 



1^2.. .3,2.. .3 1 
1^1. ..2,2.. .3 1 



|A2...3,1...2| 
|^1...2,1...2| 



Use Jacobi's Theorem and a column and row swap to obtain 

detM' = (det^) 2 " 1 ■ detM* • (-l) 1+1+3+3 
det„4 



detM' 



IAL..24...2I 
I-A2...34...2I 



detM* 
1^1. ..2,2.. .3 1 

1^2.. .3,2.. .3 1 



«2,2 



det A 



Since „4 = A\ ...3 4 ...3, we have computed the determinant of the upper 3x3 submatrix of A. This is equiva- 
lent to the element in the upper left corner of A( 2 ' in 3.1; the negatives from the cofactors in Jacobi's method 
cancel each other out. Likewise, 

• the upper right corner of A^ has the value det Ai . y.,.4; 

• the lower left corner has the value det A2...44...3; and 

• the lower right corner has the value det Ai..A2...A- 
The next (and final) condensation in Dodgson's method is 




1 1...3,1...3ll /i 2...4,2...4| 



i 2...4,1...3ll /i 1...3,2...4| 



i2...3,2...3 1 



Apply Jacobi's Theorem on A, using for M* the 2x2 interior submatrix of A, and we have 



\A\ = 



|Ai...3 / i...3||A 2 ...4 / 2...4| - |A.2...4 / l...3||Ai...3 / 2...4| 



|A.2...3,2...3| 

This is precisely the singleton element of A^>\ See Figure 3.1. 

Theorem 10. Let A be an n x n matrix. After k successful condensations, Dodgson's method produces the matrix 



( 



A+. 1+1,1. ..k+l 
^2.. 1+2,1. ..k+1 



A 



1.. 1+1,2.. 1+2 

^2.. 1+2,2.. 1+2 



.k+l,n—k...n 



\ 



/i 2..1+2,n-fc...n 
A n -k...n,n-k...n\ ) 



\ \A„-k...n,l...k+l\ \A n -k...n,2...k+2\ 

whose entries are the determinants of the {k + l) x (k + l) submatrices of A. C> 

By "successful" iterations we mean that one never encounters division by zero 
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FIGURE 3.2. Diagram illustrating how Jacobi's Theorem is applied to Dodgson's method, 
given a general n x n matrix. 




Ai+k+l,j Ai+k+l,j+k+l 



3 



Consider a (k + 1) x (fc + 1) submatrix of A, A = A'...i+it+i,/'.../'+Jt+i- Select the 2x2 minor 

M = ( Ai 'i A y+fc+i ^ . 

Its complementary minor is M* = ( A; + j ...f+jk,/+i...i+ifc ) ■ The matrix of cofactors of M in A is 



M' 



4, 



±A' 



i,j+k+l 

±A i+k+l,j A i+k+l,j+k+l 



where 



A' 



A' 



i,j+k+l 



A-i+\...i+k+\,j+\...j+k+\ 
A -i+\...i+k+\,j.,.j+k 



A' 



By Jacobi's Theorem, 



i+k+l,j+k+l 



det M' = (det^) 2-1 ■ detM* ■ (-1 



Ai...i+k,j...j+k 



or 



(3.2) 



detM 



\i+(j+k+l)+j+(i+k+l) 



4/ 




A 1 

^■i+k+lj+k+l 




A 1 

n i,j+k+\ 




A' 




Ai+i.,.i+k,j+l...j+k 





so long as the denominator is not zero. 
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Proof. We proceed by induction on k. 

Inductive Base: When k = 1, the theorem is trivial: one condensation gives 



\Ai„2,\...2\ 
I-A2...34...2I 



\M...2,2..3\ 
\M...i+2,2.-i+2\ 



\A\...l,n-\...n\ 
\ A 2..3,n-\...n\ 



\ I A„_i... tt/ i...2| I A n „i ...n,2...i+2 1 ' ' ' I A M _i. ..n,n-l...n I / 

Inductive Hypothesis: Fix A:. Assume that for all £ = 1, . . . ,/c, the ^th condensation gives us where 
for all 1 < i,j < n 



A i...i+(n-f.),j...j+(n-t) 



Inductive Step: Let i,j £ {1, ... ,n — k}. The next condensation in Dodgson's method gives us 



'J 



I - 


A i+\...i+k+\,j...j+k 


"■i...i+k,j+\...j+k+\ 


) 




A i+l...i+k,j+l...j+k 





From the inductive hypothesis, we can substitute 

/ A i...i+k,j...j+k A i+\...i+k+\,j+\...j+k+\ 

AM) 

Apply Jacobi's Theorem with 

• A = A,-,.,j_|_jt + i / y...y_|_j : _|_i 7 

• M the 2x2 minor made up of the corners of A, 

• M' the corresponding 2x2 minor of A', and 

• M* the complementary k x k minor of A 

to see that 



A 



(k+i) 



A i...i+k+\,j...j+k+\ 



(See Figure 3.2.) 



□ 



4. Why does the "double-crossing" method work? 



The new workaround is based on Theorem 10. The goal in step i of the algorithm is to compute each 
determinant A...i+fc,/.../+)c • Dodgson's method fails when the corresponding denominator of (3.2) is zero. 

However, the fraction of (3.2) is not the only way to apply Jacobi's Theorem. As long as some (i — 1) x 
(z — 1) minor of A...;+;t,/.../+)c has non-zero determinant, we can still recover, reusing most of the computa- 
tions already performed in previous steps of Dodgson's method, and calculating only a few new minors, 
again using the same approach as Dodgson's method. For example, the proof of Theorem 6 can be summa- 
rized by applying Jacobi's Theorem with 

• A of Theorem 6 standing in for A of Jacobi's Theorem; 

• the 2x2 submatrix M* of Theorem 6 standing in for M* of Jacobi's Theorem; and 

• the matrix of determinants of minors of M' of Theorem 6 standing in for M' of Jacobi's Theorem. 
One difference does require investigation: the negatives in the matrix of cofactors are now in different 
places! With Dodgson's Method, we choose the central minor for M, which has two important conse- 
quences in Jacobi's Theorem: 

• i ( = j f for each I, so that (-l) E <=i ie+ ' e = (-l)^7=i 2i = 1; and 

• if negatives appear in the 2x2 matrix of cofactors, they appear off the main diagonal, cancelling 
each other out. 

In the special case of the double-crossing method, the matter is a little more complicated. We have chosen 
the top-middle minor for M, so: 



(-1) 



Lf =1 i(+U 



_ 1 -jl+2+2+3+---+fc+fc+l _ ^ 



l) 1+fc ;and 
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FIGURE 4.1. Diagram for workaround for Example 11 





^1,3 




^3,2 




^2,4 
^3,4 



(2) (2) 

The problem in computing N^ 2 was the zero in position Na, 3 . Since N 12 



detN2...,4 / 2...4 choose 



instead the minor 



whose complementary minor is M* 



M 



/ ( _!)2+i 



M' 



(-1) 



3+1 



Nl, 3 

N 3 ,3 N 3A 
N 2 ,3 N 2/4 



N 3 ,2 N 3/4 
3. In this case the matrix of cofactors is 



2+3 



("I) 
(-1) 3+3 



/ _ 


3 






3 


\ 




2 






1 2 






3 






3 




V 


1 






-1 


/ 



Ni /2 Ni /3 

^3,2 N 3 , 3 
N 2/ 2 N 2 ,3 



By Jacobi's Theorem, 



detA/2...4 2...4 



detM' 



v 1+2+3+3 



3. 



detM* ■ (-1) J 

Notice the correspondence between M' and the matrix obtained in Example 7 using the double- 
crossing method: multiplying the top row of M' by —1 gives us N( 2 ). Dividing detM' by -3 = 
3 ■ ( — l) 9 cancels out the negative introduced into the determinant. 



the determinant of the 2x2 matrix of cofactors is 

(■_l)2t+2;+3fc-l , D 



where D is the determinant of the minors. The multiple of D simplifies to ( — l) fc 1 , which has the 
same sign as (-l) 1+,c ; thus we can disregard the signs and consider only the determinants of the 
minors. 

We use this to revisit Example 3, which would not work with Dodgson's method even after swapping rows 
or columns. 



Example 11. Recall from Example 3 



N 



/ 1 3 \ 
0-101 

1 12 

V 2 1 J 

r(2) 



Dodgson's method would have us compute AT) 2 by choosing for Jacobi's Theorem 



M 



whence 



but detM* = N 23 = 0. 



Ng =detN 1 ... 3 , 2 ...4 



M* = (N 2 , 3 ), 
detM' 



detM* 
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The double-crossing method allows us to choose instead 



M 



N 2 ,2 N 2A 
N 3 ,2 N 3A 



Since N\ 3 = 1 we have 



N- 



(2) 
1,2 



detM' 



M* = (N 1>3 ) 



detM* • (-i)^^ 3 
det M' u ■ detM' 2/2 - detM^ 2 ■ det M 2/1 
-Nl,3 

- (detM' u ■ det M' 2 2 - detM^ 2 • det M 21 ) 



(-1) 



2+1 



"Nl,3 
(_!)3+3 



3 
-1 



(-1) 



2+3 



(-1) 



3+1 



3 
1 



-9 
^3 



3. 



See Figure 4.1. 

We conclude by showing how Jacobi's Theorem likewise justifies the "double-crossing" method in Ex- 
ample 5. 

Example 12. Recall from Example 5 



/ 1 

1 


V 2 



1 \ 

1 


1 J 



When computing Dodgson's method wants us to cross out the middle three rows and the middle 

three columns of A, obtaining the minor 



M 



n,5 



A 5 ,i A 5/5 



A 2 ,2 


^2,3 


A 2A \ 






3 


1 


^3,2 


^3,3 


Ma 


1- 




2 


1 


A 4 ,2 


^4,3 








1 


1 



whose complementary minor is 



M* 



Unfortunately the determinant of this latter matrix is zero, which corresponds to the zero in the interior of 

/ -4 -2 2 
A^ =1-2 0-2 
V 4 -4 3 

The double-crossing method suggests instead to cross out the top three rows and middle three columns 
of A, obtaining the minor 

^4,1 A 4 , 5 \ = (0 
A 5/1 A 5>5 I I 2 1 



M 



whose complementary minor is 



A u A 1/3 A 1/4 
M* = I A 2>2 A 2>3 A 2A 
A 3 ,2 A 3 ,3 A 3A 







1 




)- 


i 


3 








2 
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The determinant of M* is A({ = -2; we will multiply it by (_i) i +^+^+^+^ ) + 4 = _\ j com pute A\J we 
need det M', which requires us to compute the matrix of cof actors of M in : 



M' = 



/ 
/ 


U 


-1 
1 


u 


1 
1 




1 
1 


n 
U 


-1 
1 


U 






5 


3 


1 










5 


3 


1 






3 


2 


1 


1 




1 


3 


2 


1 









2 





1 




2 





2 












1 





1 




1 





1 









5 


3 


1 










5 


3 


1 






3 


2 


1 


1 




1 


3 


2 


1 




V 


1 


1 


1 










1 


1 


1 


/ 



The top row of negatives introduces a negative into det M', which cancels with the — 1 that is multiplied to 
det M* , cancelling each other out. We can thus consider the determinants of minors, rather than the matrix 
of cof actors. O 

Of course, it is possible that both an element of the interior is zero and the element above it is zero. The 
following theorem provides the promised generalization of Theorem 6; its proof uses Jacobi's method in 
a manner similar to the proof of Theorem 6. Notice that if r — s — 0, then Theorem 13 specializes to 
Dodgson's method. 

Theorem 13. Let A be an n x n matrix. Suppose that we try to evaluate \A\ using Dodgson's method, hut we 
encounter a zero in the interior of A^ k \ say in row i and column j of A^ k \ Let r, s e { — 1,0, 1}. If the element a in 
row i + r and column j + s of is non-zero, then we compute A^ 2 ^ as usual, with the following exception for 
the element in row i — 1 and column j — 1: 

• let £ = n — k; 

• identify the {£ + 3) x (£ + 3) submatrix A whose upper left corner is the element in row i — £ and column 
j-icfAW; 

• identify the complementary minor M* by crossing out the (I + 1) x (I + 1) submatrix M of A whose upper 
left corner is the element in row r + 2 and column s + 2 of A; 

• compute the matrix M' of determinants of minors ofM* in A; 

• compute the element in row i — 1 and column j — 1 of A( fc+2 ) by dividing the determinant of M' by a. 
We can use Dodgson's method to compute the determinants ofM'. Q> 

5. Concluding remarks 

The double-crossing method succeeds as long as Dodgson's method succeeds, and with the same num- 
ber of integer operations, since Dodgson's method is a special case of the double-crossing method. In the 
worst case scenario where Dodgson's method generates a number of interior zeroes, the double-crossing 
method requires the computation of at least two (k — 1) x (k — 1) determinants for every zero entry of A^ k \ 
However, this changes the number of integer operations only by a constant, so the double-crossing method 
requires only O (n 3 ) integer operations. In addition, the double-crossing method preserves the general 
simplicity and spirit of Dodgson's method. 

The double-crossing method is not guaranteed to succeed; if an intermediate matrix contains a 3 x 3 
block of zeroes, then the double-crossing method also fails. Sparse matrices provide an excellent example 
where the double-crossing method is an abject failure; consider the identity matrix of order six or higher. 

In those cases where the double-crossing method fails, one can still preserve the computations that work, 
and adapt a hybrid with another method. Theorem 13 tells us the precise minor A whose determinant we 
need; we can compute the determinant of this minor using another method, substitute its value into row 
i — 1 and column; — 1 of A^ k+2 \ and proceed. 
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